Accelerated Galerkin BEM for 3D Potential Theory
نویسنده
چکیده
This paper is concerned with the development of a fast spectral method for solving boundary integral equations in three-dimensional potential theory. Upon discretizing the underlying boundary integrals via a Galerkin approximation, the proposed method overlays the problem domain with a regular Cartesian grid that serves as an auxiliary platform for computation. With the aid of the Fast Fourier Transform, the necessary influence matrices of the discretized problem are rapidly evaluated on the regular grid in a sparse manner. Unlike traditional techniques dealing with boundary integrals, the sparse representation of the featured coefficient matrices results in a significant reduction in computer memory requirements. The computational cost associated with the sparse approximation of influence matrices is asymptotically lower than that of conventional methods. For a numerical solution of the resulting linear system, a Krylov-subspace (e.g. BiCGSTAB) iterative method is further employed wherein the sparse influences are used to rapidly compute the matrix-vector products involved at each iteration. Several key features of the formulation, including the mapping of density distributions onto the regular Cartesian grid, are highlighted. Numerical experiments are presented to illustrate the performance of the spectral method. The proposed approach will find application in areas involving large simulations with complex and moving boundaries.
منابع مشابه
A pFFT Accelerated Linear Strength BEM Potential Solver
A linear strength, Galerkin Boundary Element Method (BEM) for the solution of the three dimensional, direct potential boundary integral equation is presented. The method incorporates node based linear shape functions of the single and double layers on °at triangular elements. The BEM solution is accelerated using a precorrected Fast Fourier Transform algorithm (pFFT) [1]. Due to the extended co...
متن کاملExploiting partial or complete geometrical symmetry in 3D symmetric Galerkin indirect BEM formulations
متن کامل
Classical FEM - BEM coupling methods : nonlinearities , well - posedness , and adaptivity
We consider a (possibly) nonlinear interface problem in 2D and 3D, which is solved by use of various adaptive FEM-BEM coupling strategies, namely the Johnson-Nédélec coupling, the Bielak-MacCamy coupling, and Costabel’s symmetric coupling. We provide a framework to prove that the continuous as well as the discrete Galerkin solutions of these coupling methods additionally solve an appropriate op...
متن کاملAdaptive Multilevel BEM for Acoustic Scattering
We study the hand p{versions of the Galerkin boundary element method for integral equations of the rst kind in 2D and 3D which result from the scattering of time harmonic acoustic waves at hard or soft scatterers. We derive an abstract a{posteriori error estimate for indeenite problems which is based on stable multilevel decompositions of our test and trial spaces. The Galerkin error is estimat...
متن کاملCoupling of Raviart-Thomas and Hybridizable Discontinuous Galerkin Methods with BEM
We present new a priori error analyses of the coupling of the Raviart–Thomas (RT) method and BEM as well as of the coupling of the hybridizable discontinuous Galerkin (HDG) method with BEM. The novel features of the analysis of the RT-BEM coupling are the superconvergence estimates for the scalar approximation of the interior problem and new bounds in weak and strong norms for the boundary vari...
متن کامل